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Quantum Heisenberg antiferromagnets on pyrochlore and checkerboard lattices in a strong exter- 
nal magnetic field are mapped onto hard-core lattice gases with an extended exclusion region. The 
effective models are studied by the exchange Monte Carlo simulations and by the transfer matrix 
method. The transition point and the critical exponents are obtained numerically for a square-lattice 
gas of particles with the second- neighbor exclusion, which describes a checkerboard antiferromagnet. 
The exact structure of the magnon crystal state is determined for a pyrochlore antiferromagnet. 

PACS numbers: 75.50.Ee, 75.45.+j, 75.10.Jm, 75.30.Sg 



I. INTRODUCTION 



A characteristic feature of geometrically frustrated 
spin models is an infinite degeneracy of the classical 
ground state. Such a degeneracy may be lifted by quan- 
tum or thermal fluctuations via the 'order from disor- 
der' effect jii^i^ A closely related property is emergence of 
small energy scales, which are described by the effective 
Hamiltonians acting in the subspace of nearly degener- 
ate low-energy quantum states. To date, several types of 
effective models have been developed for frustrated an- 
tiferromagnets in different regimes. The first example 
is a celebrated quantum dimer model introduced phe- 
nomenologically by Rokhsar and Kivelson to describe a 
magnetically disordered resontaing valence bond phase of 
S = \/2 antiferromagnets4 Recently, effective quantum 
dimer models have been derived for several realistic spin 
Hamiltonians in zero and in a finite magnetic fieldi^i^ii 
The second example is the effective Hamiltonian for the 
collinear states, which applies to a semiclassical {S ^ 1) 
pyrochlore antiferromagnet in zero magnetic field.™ The 
third type of effective models describes geometrically 
frustrated magnets at high magnetic fieldsi^ '^°i^^'^^ These 
are the so called lattice-gas models considered further in 
the present work. 

Previously, the effective lattice-gas models have been 
derived and discussed for the sawtooth chain and a 
kagome ant iferromagneti ^'^^i^^ These two quantum an- 
tiferromagnets are mapped onto lattice gases of hard- 
core classical particles with the nearest-neighbor exclu- 
sion on a one-dimensional chain and a triangular lat- 
tice, which are respectively called the hard-dimer and the 
hard-hexagon model. In both cases, the thermodynamics 
of the effective models has been calculated exactlyJ^ In 
the present article, we shall focus on the high-field be- 
havior of quantum antiferromagnets on two- and three- 
dimensional pyrochlore lattices. The corresponding effec- 
tive models do not allow exact solution for their thermo- 
dynamic properties, therefore, we shall study them via a 





(a) 



(b) 



FIG. 1: Pyrochlore lattice of corner-sharing tetrahedra (a) 
and its two-dimensional analog, checkerboard lattice (b). 
Solid lines show the nearest-neighbor bonds of strength J. 
Dashed lines indicate possible further neighbor exchanges. 



combination of numerical techniques. 

The pyrochlore and the checkerboard lattices (Fig. [T|) 
consist of 3D and 2D networks of corner-sharing tetrahe- 
dra (four-spin blocks). For the most part of the present 
work, we disregard weak further neighbor exchanges and 
consider the nearest-neighbor Heisenberg model with an 
arbitrary spin S in an external magnetic field: 



H 
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Above the saturation field Hg the ground state of a 
Heisenberg antiferromagnet is a fully polarized vacuum 
|0) = I ITT ■■■), where all spins are in a state with the 
maximum possible value of Sf = S. The low-lying ex- 
citations are single spin fiips \i) = S^\0), whose disper- 
sion can be straightforwardly calculated. The charac- 
teristic feature of many geometrically frustrated models 
including the two considered models is presence of the 
low-energy branch of localized magnons with the energy 
eo = The details for the two discussed 
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models are presented in the subsequent sections. 

Localized nature of magnons from the lowest branch 
allows to construct a subclass of exact multimagnon 
statesjiS a quantum state with m localized magnons (m- 
LMs) occupying spatially separated regions of the lattice, 
which are not directly connected by the exchange bonds, 
is an exact eigenstate of the quantum Hamiltonian ([T]) 
and its energy is given by meo. It has been proven ex- 
actly that isolated LM states correspond to the lowest 
energy states in every magnetization subsectorj -'^^i^^ In 
addition, the renormalization arguments suggest that de- 
localized (propagating) states are separated by a finite 
gap from LMs»2.iiS Hence, at low temperatures the parti- 
tion function of a quantum frustrated antiferromagnet in 
the vicinity of the saturation field is determined entirely 
by LMs: 

J2 ^(™' ^) ^"^^'^ ' n = Hs-H , (2) 

m=0 

with z = e'^l'^ being the activity. Combinatorial factors 
wirn^ N) denote a number of linearly independent m-LM 
states, while A'max is the maximal number of LMs for 
a lattice with N sites {N ^ oo). Mapping to a lattice 
gas of hard-core objects is used as an approximate way 
to calculate z«(m, N) and Z. The general consequence 
of Eq. ^ is the following scaling of the total entropy: 
S = f[{H - Hs)/T]. At H = Hs the entropy is temper- 
ature independent and, consequently, even at T = a 
frustrated quantum antiferromagnet has a finite macro- 
scopic entropy, which is determined purely by the lattice 
geometry. The quantum order from disorder mechanism 
becomes ineffective for this special value of applied field 
due to the localized nature of (exact) quantum states 
with zero energy. In Sees. II and III we describe deriva- 
tion of lattice gas mapping and obtain conclusions sep- 
arately for the checkerboard and the pyrochlore antifer- 
romagnet. Role of extra exchanges beyond the nearest- 
neighbor pairs, see Fig. [2 is briefly discussed in Sec. IV. 

II. CHECKERBOARD ANTIFERROMAGNET 
A. Localized magnon states 

A checkerboard lattice contains two spins in the prim- 
itive unit cell. The corresponding Bravais lattice is 
formed by elementary translations on ai — (1,1) and 
a2 = (—1,1) and has a square shape. In the saturated 
phase at high fields, the one-magnon spectrum has, ac- 
cordingly, two branches: 

t^ik = - SJS , uJ2^. = H- US{1 - 7k) , (3) 

where 7k = cos cos ky . The saturation field Hs — SJS 
corresponds to the vanishing energy of magnons from the 
lowest dispersionless branch wik- Localization of excita- 
tions from the lowest branch is determined by the lattice 
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FIG. 2: (color online). Localized magnons on a checker- 
board lattice. Thick lines indicate positions of one-magnon 
states with corresponding phases shown by -I- and — . The 
shaded area shows excluded volume around a central square- 
void state, which has to be respected for construction of mul- 
tiparticle states. Four LMs in the upper-left part satisfy a 
close-packed condition (see the text). 

topology. A simple localized state can be constructed as 
a spin-flip trapped on a square void of a checkerboard 
lattice: 

|^.) = -L^(-l)"-i^-,|0), (4) 

^ n— 1 

where the numbering of sites goes counterclockwise start- 
ing from the lowest right corner, see Fig. [5) The prob- 
ability to find spin-flip on a site adjacent to the void 
vanishes due to the destructive interference. This prop- 
erty crucially depends on equal strength of all bonds of 
the checkerboard lattice. 

A linear combination of two square-void states shar- 
ing one vertex encompasses two square voids and so on, 
see Fig. [21 An arbitrary 1-LM state can be constructed 
by drawing a closed graph on the original lattice, which 
passes through two vertices of every crossed four-site 
block, and assigning -I- and — signs in alternate order for 
a spin-flip propagating around it. The smallest square- 
void states play an important role by forming a com- 
plete nonorthogonal basis in the subspace of dispersion- 
less one-magnon states^iS, The multi-particle LM states 
are, then, constructed by respecting topology of the ex- 
change bonds: LMs cannot occupy square voids, which 
are contiguous to the same crossed square as shown in 
Fig. [2] by shaded area. 

There is an upper limit on the density of isolated local- 
ized magnons. For a checkerboard lattice, a close-packed 
structure is constructed by putting LMs on every sec- 
ond square void in horizontal rows with an alternating 
shift between the rows, see the left part in Fig. O This 
pattern corresponds to a magnon crystal, which breaks 
translational symmetry and has the density of LMs equal 



to no = iVmax/-^ = 1/8. The breaking of the transla- 
tional symmetry is, however, incomplete: diagonal rows 
of LMs can freely slide without affecting magnons in ad- 
jacent rows. The degeneracy of the magnon crystal state 
is A^dog ~ 2^+^, where i is a linear size of the system, see 
also Sec. lie. Note, that the magnon crystal for a quan- 
tum kagome antiferromagnet is only three-fold degener- 
ate. In that case there is a well-defined finite tempera- 
ture transition associated with the translational symme- 
try breaking. 

The purpose of the present study is to investigate the 
low-temperature behavior of the checkerboard antiferro- 
magnet and nature of a phase transition into the magnon 
crystal state. In order to proceed, we map a quantum 
checkerboard antiferromagnet in the vicinity of the sat- 
uration onto a gas of hard-core particles. Such an effec- 
tive model is defined on a square lattice formed 
by centers of square voids of the original checkerboard 
lattice. The dual square lattice is rotated by 45° and 
contains N/2 sites. Localized magnons are represented 
by hard-core classical particles obeying the nearest- and 
the next-nearest neighbor exclusion principle. Below, we 
call them 'hard-polygon states.' One such polygon for 
a checkerboard lattice is shown in Fig. [2] by a shaded 
square. In reality a state with two localized magnons oc- 
cupying adjacent sites is separated by only a finite gap 
from the low-lying LM states. An estimate of such gap 
for a kagome lattice antiferromagnet is given in Ref.H In 
the following only T — > regime is considered, where con- 
tribution of such higher energy states can be neglected. 



B. Topological classes of localized magnons 

The hard-polygon states being linearly independent^^ 
do not exhaust all possible localized magnon states. 
The missing 'defect' states belong to different topolog- 
ical classes of LMs. 2'— Let us consider an infinite plane 
with open boundary conditions, which has topology of 
a sphere without the North pole. As was explained be- 
fore, one-particle localized states correspond to closed 
lines. Successive expansion of the wave- function of a 
one-magnon state in terms of basis plaquette states can 
be represented as a gradual deformation of a long loop 
with subsequent contraction into a point (center of the 
last plaquette). Two closed loops for a 2-LM state may 
be either contractible to two distinct points on a plane 
or lie inside each other and, hence, be contractible into 
a single point, see Figs. [3^ and [31d. In the first case 
the two magnon state belongs to a subset of two-particle 
hard-polygon states, whereas in the second case an ex- 
pansion in the hard-polygon states is impossible. The 
two-magnon states shown in Figs. [3^ and [3)d are, there- 
fore, linearly independent. 

Situation is somewhat different for a cluster with peri- 
odic boundary conditions, which has a torus topology. In 
that case one can continuously deform the outer line in 
Fig.Oa, split it into two loops with finite winding around 
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FIG. 3: (color online). Examples of topologically different 
states of localized magnons: 2-magnon subsector (a) and (b); 
3-magnon subsector (c)-(e); (f) shows an elementary defect 
state from the class (b). 



the torus, move them around and close again into one 
contour on the opposite side such that the two-magnon 
state in Fig. transforms into a state of type Fig. [3^. 
The linear dependence of two-magnon states in Figs. [3^ 
and [3)3 for the torus topology is explained by presence in 
this case of one linear relation between 1-LM square- void 
(hexagon) states Still, there are additional 2-LM states 
in this topology, which are constructed by putting one 
or two LMs on contours with nontrivial winding around 
the torus. Topologically nontrivial classes of LM states 
in the three-magnon subsector for a lattice with open 
boundary conditions are shown in Figs. [3J;-e. For peri- 
odic boundary conditions the state in Fig. ^ becomes 
topologically equivalent to the state in Fig. [5^, while the 
state Fig. [3jl can be transformed into the hard-polygon 
state. The macroscopic limit, however, does not depend 
on the boundary conditions. Equivalence of the two ap- 
proaches in the — > 00 limit is recovered by observing 
that the largest number of topologically nontrivial LM 
states in n-magnon subsectors with n > 2 is given by 
states of the type shown in Fig. [3j3, which are present for 
both choices of boundary conditions. Disentanglement of 
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two, three, etc. enclosed loops in the torus topology is im- 
possible once other LMs are present. On the other hand, 
the contribution from closed loops with a nontrivial wind- 
ing around the torus corresponds in the thermodynamic 
limit only to a surface effect. 

Topological origin of the additional localized magnon 
states determines their presence for all two-dimensional 
frustrated models including antiferromagnets on kagome, 
checkerboard, and star lattices. There are no such states 
in one-dimensional models as, for example, the sawtooth 
chain, where the hard-particle representation is asymp- 
totically exact. In order to estimate the contribution 
of additional LMs, one has to define the basis states in 
the topologically nontrivial classes. The elementary de- 
fect state in the topological class of Fig. [Sb is shown 
in Fig. [3f (open boundary conditions are assumed). Its 
wave-function is given by 

8 

|2-defect) c V«> + l^o>) l^o) , (5) 

4 = 1 

where the numbering of plaquette states follows Fig. [3f. 
This wave- function includes two LMs on adjacent voids 
and on the central plaquette and violates, therefore, the 
hard-polygon constraint. The two-magnon states ([5]) re- 
siding on different central squares are linearly indepen- 
dent and form a nonorthogonal basis in the subspace of 
topologically nontrivial graphs of Fig. [3}d. The 2-LM de- 
fect state ([U can be identified with a new classical par- 
ticle, which has energy 2eo and a longer-range repulsion. 
Such a mapping suggests that the topologically nontriv- 
ial states yield an additional macroscopic contribution to 
the partition function Z. Their share is, however, sup- 
pressed compared to the hard-polygon contribution by a 
large entropic factor: region occupied by the basis two- 
magnon state. Fig. [3f, can be occupied by 16 different 2- 
LM states residing on small empty squares. In addition, 
by counting the total number of excluded square voids 
we conclude that the states (O do not contribute appre- 
ciably to the magnetization subsectors with the average 
density of LMs larger than n > 2/(25 • 2) = 0.04.^^ Con- 
sequently, their role at high densities, e.g.^ in the vicinity 
of the transition into a magnon crystal state is strongly 
suppressed. Alternatively, the lattice gas mapping can be 
improved by including extra types of particles, which de- 
scribe LMs from different topological classes. Below we 
consider the simplest version of the lattice gas mapping 
with only one type of particles. 

C. Effective lattice-gas model 

Two-dimensional lattice gas models were originally 
suggested to study atomic adsorption on various 
substrates.^^'22^'2i^'26^^ They also describe low- 
T behavior of Ising antiferromagnets in a longitudi- 
nal field] ^^'^^1'^° Despite the long-lasting interest there 
are only few well-established facts for a square lattice 



TABLE L The entropy per site <S and the density n of a 
hard-square lattice gas with second-neighbor exclusion at /i = 
obtained from the transfer matrix calculation on a semi- 
infinite strip A'f X oo. 
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0.13713 
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0.294647 


0.13679 


14 


0.294642 
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gas with the nearest- and the next-nearest-neighbor ex- 
clusion. It is generally accepted that this model ex- 
hibits a continuous phase transition into a partially or- 
dered (2x1) phase This phase transition belongs 
to the universality class of an XY model with a four- 
fold anisotropy,-- which is a special case of a more gen- 
eral Z{A) discrete planar model.— i^i^ Depending on the 
values of two coupling constants the Z{A) model exhibits 
either two Ising-like transitions between a paramagnetic 
and an ordered state or a single critical point, "^^ which has 
nonuniversal critical exponents.'^'* Estimates for the tran- 
sition point of a hard-square lattice gas with the second- 
neighbor repulsion vary from Inz^ = (/Lt/T)c = 5.3 in the 
early work^^ to In Zc = 4.7 in the later studyi^ There is 
even less certainty about values of the critical exponents. 
The entropy of the system at /i = 0, which is an interest- 
ing property of a quantum antiferromagnet, has not been 
investigated in the context of previous applications. 

First, we investigate numerically the entropy of the lat- 
tice gas model at zero chemical potential. This can be 
done by the transfer matrix method and Monte Carlo 
simulations. The former method has been a standard 
technique since the early studies of two-dimensional lat- 
tice gas models^ii^^i^ In this scheme a two dimensional 
lattice is represented by a semi-infinite strip of width M 
and the thermodynamic properties are derived from the 
largest eigenvalue Ai of the transfer matrix. In partic- 
ular, the entropy normalized per one site is equal to 
S = (lnAi)/Af. The actual calculations become quite 
simple for /i = 0, when the elements of the transfer ma- 
trix take only values or 1 . The results for a few values 
of M are presented in table I. The convergence is very 
rapid and already the M — 10 strip gives four significant 
digits for the entropy. 

To determine entropy from Monte Carlo (MC) simu- 
lations we adopt the following procedure. The standard 
Metropolis algorithm is used for gradual annealing from 
the low-activity regime Inz — —20, where the density of 
particles is vanishingly small and 5 = 0, to the point 
Inz = 0. The step for the chemical potential is chosen 
to be small enough A(ln2;) = 0.05. At every value of fi 
10^ MC steps (lattice sweeps) are performed for equili- 
bration and after that 10^ MC steps are used to measure 
{dS / d , which is calculated from the cumulant of the 
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FIG. 4: (color online). Effective lattice gas model for a 
checkerboard antiferromagnet. The dual square lattice is 45° 
rotated compared to the original checkerboard lattice in Fig. 2 
and has twice less sites. Hard particles are shown by large 
light circles; small dark circles indicate the excluded sites 
around one particle. The presented close-packed structure 
with half-filled rows randomly shifted in the horizontal direc- 
tion has a nonzero Fourier harmonic at Qi = (0, tt). 

energy E and the number of particles M: 

In our case E — ~^Af and the cumulant in Eq. ^ is 
proportional to the variance of the total number of par- 
ticles. The error bars are estimated by performing up 
to 100 independent runs with different sequences of ran- 
dom numbers. Afterwords, the data for {dS/d^)T are 
numerically integrated to find iS|^=o- Simulations have 
been performed on square clusters with periodic bound- 
ary conditions and N = sites, L = 16, 32, and 64. The 
obtained results agree with each other within numerical 
accuracy and yield S = 0.2946(1) for the entropy per site 
and n = 0.13676(1) for the average density of particles. 
The found values are in good correspondence with the 
transfer-matrix results included in table I confirming the 
accuracy of both methods. The particle density at = 
is quite substantial and only two times smaller than the 
density of the close-packed structure ng = 0.25. 

Once the chemical potential further increases (z > 1) 
more and more particles become condensed. Eventually 
the lattice gas transforms into an ordered close-packed 
structure. The order parameter of such a crystalline 
state are certain Fourier harmonic(s) of the particle den- 
sity, riq = (1/A^) n^e"*'!'"'. A close-packed structure 
with half- filled rows, which are randomly shifted in the 
x-direction has a nonzero harmonic at Qi = (0,7r) with 
riQj = 1/4, see Fig. |4l A similar state with random 
shifts of columns along the y direction is described by 
Q2 = (tt, 0) and hq^ = 1/4. These two phases are called 
(2 X 1) states in a standard nomenclature adopted for lat- 
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FIG. 5: (color online). Dependence of the two order parame- 
ters Eq. (O on logarithm of the fugacity In z = fi/T for three 
system sizes. 

tice gas models.™ A fully symmetric structure with all 
rows (columns) in phase with each other has equal ampli- 
tudes uq. = 1/4 for the three wave- vectors Qi, Q2, and 
Q3 = (tt, tt) and is called a (2x2) state. The wave- vectors 
Qi and Q2 belong to the same irreducible representation. 
Consequently, we consider two order parameters: 

^l = ("Q, +"qJ'^' ' M2 = M. (7) 

The infinitely degenerate close-packed structure for 
a square-lattice gas with the second-neighbor exclusion 
leads to equilibration problems in simple MC simulations. 
This was evident since the early MC study of a frustrated 
Ising antiferromagnetf^l which investigated clusters up to 
only 40^ sites using a single spin-flip relaxation method. 
We employ instead an exchange MC algorithm proposed 
some time ago to tackle systems with extremely long re- 
laxation times, as, for example, spin glassesi^^ With this 
modification we have been able to study clusters with up 
to 120^ sites. To ensure a substantial replica exchange 
rate ~ 0.7-0.9, we have used 120 replicas in the range 
< Inz < 8 for all system sizes. The simulation runs in- 
cluded 10^ exchange MC steps for equilibration and up to 
10^ MC steps for measurements. Statistical errors were 
estimated from bining the MC series for each value of z. 

The specific heat exhibits a broad and rounded max- 
imum near to In z ~ 4.4. The height of the peak grows 
very slowly with increasing L, such that a preliminary de- 
termination of the transition point is impossible from the 
specific heat data. The ensemble averages of the two or- 
der parameters (Mi) and (M2) are shown in Fig. [51 The 
square of the second order parameter (M|) goes down to 
zero as 1/A^ with the system size, whereas (Mi) scales 
to a finite constant at Inz > 5. The limiting value for 
the largest cluster is very close 1/4 in accordance with 
the previous analysis for the (2x1) type of ordering. 
The precise location of the transition point is obtained 
by measuring the Binder cumulant U4 — {M^) / {Ml)"^ . 
Its dependence on z is presented in Fig. [5] The tran- 
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Inz 

FIG. 6: (color online). Binder cumulant U4 as a function 
of Inz = fi/T for different cluster sizes L. The inset shows 
scaling of the derivative dU4/d{lnz) at the transition point 
In Zc — 4.56 versus L. 

sition point can be estimated from the crossing points 
of Binder cumulants for different clusters. '^^ This yields 
Inzc — 4.56 ± 0.02 for the critical activity. The density 
of particles at the transition point is Uc = 0.2325(3). 

In principle, the measurement of Mi ^ cannot dis- 
criminate between a single-fc structure corresponding to 
the (2x1) order and a double-fc structure with two 
nonzero Fourier harmonics nq^ , nq^ ^ 0. The two struc- 
tures can be distinguished with the help of an order pa- 
rameter, which probes breaking of the C4 rotational sym- 
metry: 

A = ^ n,{n.,+2^ - n,+2y) ■ (8) 

i 

The order parameter A vanishes for the (2 x 2) state and 
all other states symmetric under 90° rotations, whereas 
it has a finite value in the (2x1) phase. We have cal- 
culated (|A(2;)|) by MC simulations and found that at 
In z > 5 within statistical errors it reaches a finite value: 
(|A|) = 1/8. This not only demonstrates an absence of 
the tetragonal symmetry in the ordered state, but also 
proves random shifts of half- filled rows (columns). For 
example, if shifts occur in a regular alternating order, 
than one would find (|A|) — 1/4, which is definitely ex- 
cluded by our Monte Carlo results. 

In order to determine the critical exponents, we per- 
formed finite-size scaling analysis of various thermody- 
namic quantities measured at the estimated critical point 
In Zc- The correlation length exponent v is extracted from 
the behavior of the derivative of the Binder cumulant 
dUi/d{lnz) ^ L^/'^ at z = Zc- The fitting shown in the 
inset of Fig. [H yields l/v = 1.16 ± 0.02. The order pa- 
rameter at criticality scales as M{zc) ^ L~^/'^ . A fit for 
a few largest L gives /3/^ = 0.15(1). As a result, we ob- 
tain the following estimates v = 0.86(2) and /3 — 0.13(1). 
The largest clusters employed in the present study are 



still not sufficient to independently extract the critical 
exponent a, presumably because of a large regular con- 
tribution to the specific heat compared to a universal 
singular part. The obtained values for /3 and v place a 
hard-square gas with the second-neighbor exclusion on 
a line of critical points of the Z{A) model between the 
Potts model, which has /3 = 1/12 and v = 2/3,— and the 
vector Potts model, which belongs to the Ising universal- 
ity class with /3 = 1/8 and v = l.'^^ Our result for Inzj, 
coincides within the error bars with the corresponding 
value from the recent independent MC study of the same 
model^- though we do not share its claim of the Ising 
universality class for the transition. 

Properties of a quantum spin-1/2 checkerboard anti- 
ferromagnet in the vicinity of the saturation field are ob- 
tained from the above results by a straightforward rescal- 
ing to the number of sites of checkerboard lattice, which 
is twice larger than the number of sites for the lattice gas. 
Ax H ~ Hg, the entropy and the magnetization normal- 
ized per one site have universal T-independent values: 

Ss = 0.1473 , M, = 0.4316 . (9) 

Note, that the entropy of a checkerboard antiferromagnet 
is 30% larger than the corresponding result for a kagome 
antiferromagnet.^ The transition field into the magnon 
crystal phase is given by 

iJc(T) = iJ, -Tlnzc = 4J-4.56T . (10) 

Let us emphasize again that the entropy found within 
the lattice-gas (hard-polygon) description of LMs is only 
approximate. Localized magnons from different topolog- 
ical classes will increase the above value of Sg- Such 
corrections should not be very large in view of a large 
size of basis defect states though it would be in- 
teresting to estimate ASs from an improved lattice gas 
mapping, as discussed in the end of Sec. lib, or from 
the exact-diagonalization data. Propagating n-magnon 
states can be separated by rather small gaps from n-LM 
states. Therefore, Eq. PH)) describes the slope of the 
actual transition line Hc{T) in the frustrated quantum 
antiferromagnet at T ^ 0. 

III. PYROCHLORE ANTIFERROMAGNET 
A. Localized magnons 

The primitive unit cell of a pyrochlore lattice contains 
4 spins belonging to one tetrahedron. These tetrahe- 
dra are arranged into an fee lattice, which is formed 
by the elementary translations on ai — (0,1/2,1/2), 
&2 = (1/2,0,1/2), and ag = (1/2,1/2,0). (Linear size 
of the standard cubic cell with 16 spins is chosen as the 
unit of lengthi^^) In a strong magnetic field, where the 
saturated phase is stable, the four one-magnon excitation 
branches are given by 

2 = H- 8JS, UJ3A = H~ US ± 2JSy/T+^, 
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FIG. 7: (color online). Cross-section of a pyrochlore lattice 
perpendicular to the [111] axis. Tetrahedra from the same 
kagome layer are drawn by solid lines. Thick line shows a 
localized magnon on one of the hexagons in the kagome plane. 
Dashed lines indicate tetrahedra from the upper kagome layer 
(shown only partially). Dots with numbers denote centers of 
hexagons in horizontal and tilted planes, which form a dual 
pyrochlore lattice. 

rju = cos cos — + cos — cos h COS COS . 11 1 

' 2 2 2 2 2 2 ^^ 

The saturation field of an antiferromagnet on a py- 
rochlore lattice is the same as for a checkerboard lattice: 
Hs ~ 8JS. Two dispersionless branches 0^1,2 contain 
2 X (A^/4 - 1) = N/2 - 2 linearly independent local- 
ized modes, N being the number of pyrochlore lattice 
sites. Geometric interpretation of these modes is essen- 
tially the same as for two-dimensional frustrated lattices. 
The smallest LMs are located on hexagon voids of kagome 
layers, which alternate with triangular layers along the 
[111] and three other equivalent directions, see Fig. [T] 
An A''-site pyrochlore cluster contains N hexagons, which 
for A'^ — > 00 is twice more than the number of localized 
modes in the two lowest branches ([TTT) . The 'elementary' 
hexagon modes are, therefore, not only nonorthogonal 
but also form an overcomplete basis. This can be ex- 
plained by observing that hexagon loops on a pyrochlore 
lattice obey exactly N/2 linear relations, which leaves 
only N/2 linearly independent stateSi^- Despite this fact 
we shall use below the hexagon states as a basis in a map- 
ping of a pyrochlore antiferromagnet near on a lat- 
tice gas model. Such a representation gives an incorrect 
(larger) number of states in the magnetization subsectors 
with nonmacroscopic numbers of LMs. The mapping be- 
comes not so bad for the multi-magnon subsectors in the 
vicinity of close-packed structures, when the allowed con- 
figurations do not include linearly dependent states. In 
particular, using the lattice gas representation we shall be 
able to describe an exact structure of the magnon crystal 
for a pyrochlore antiferromagnet. 

The lowest-energy multi-particle states are constructed 
using the same prescription as for two-dimensional frus- 
trated models: LMs are represented by closed even-sites 



graphs, which are not directly connected by the exchange 
bonds. For example, if a LM occupies a hexagon void 
denoted by 1 in Fig. [3 then no LMs can be placed on 
hexagons denoted by numbers 2, 3 and 4, since they ei- 
ther share sites with the hexagon 1 or are connected to it 
by the nearest-neighbor bonds. The lattice gas model 
is constructed on a dual lattice formed by centers of 
hexagons of the original pyrochlore lattice, such that an 
exclusion principle for neighboring sites reproduces the 
above rules for LMs. Remarkably, the dual lattice is 
again a pyrochlore lattice. (This explains, in particu- 
lar, why the number of hexagons is equal to the number 
of sites.) If the nearest-neighbor distance is denoted by 
d {d=V2/4), then presence of a particle on a given site 
excludes for occupation (i) 6 nearest-neighbor sites at the 
distance d as, for example, a pair of sites (hexagons) 2 and 
4 in Fig. [3 (ii) 12 second-neighbor sites at the distance 
\/3d, e.g., a pair of sites 1 and 2, (iii) 6 third-neighbor 
sites at the distance 2d, e.g., sites 1 and 3, and (iv) 12 
fourth-neighbor sites at the distance ^/Ed, e.g., sites 1 
and 4. Note, that on a pyrochlore lattice there are two 
types of third-neighbor sites at the distance 2d. The first 
type used above corresponds to pairs on opposite vertices 
of hexagons. The second type of third neighbors corre- 
spond to pairs with a third site in between. There are 
again 6 third-neighbor sites of this type around a given 
site. In terms of the original lattice they correspond to 
hexagon pairs (2,5) and (2,6) in Fig. [71 which allow si- 
multaneous occupation by LMs. These hexagons belong 
to two different parallel kagome planes. 

Let us now consider possible magnon crystal states ob- 
tained by a close-packing of the smallest LMs on hexagon 
voids. The initial estimate for the maximal density of 
LMs has been obtained from the known result for a sin- 
gle kagome planeJ^ filling one-third of hexagons in every 
kagome plane, which is perpendicular to one of the four 
cubic diagonals, gives uq — (1/3) x (1/4) — 1/12 for the 
density. No arguments were given yet that this is indeed 
the maximal possible density. In terms of the effective 
lattice gas model, which operates on a dual pyrochlore 
lattice, a close-packed structure is constructed by fill- 
ing every second site in one of the six chains formed 
by six edges of one tetrahedron, see Fig. [8l Particles 
at the distance 2d along the chains correspond to the 
second type of third-neighbors and are allowed by the 
exclusion principle. All chains parallel to a chosen direc- 
tion form a triangular lattice in the perpendicular plane. 
Half-filling of the nearest-neighbor chains is prohibited by 
the exclusion principle, whereas second-neighbor chains 
in the triangular plane allow for simultaneous occupa- 
tion. In this way it is possible to put particles only 
on one-third of the parallel chains, which again yields 
no = (1/2) X (1/3) X (1/2) = 1/12 for the particle density 
(The last factor 1/2 corresponds to the fact that parallel 
chains contain only a half of all lattice sites, the other 
half belongs to the perpendicular chains.) An example 
of the constructed structure is shown in Fig. [51 Particles 
in every chain can be independently shifted by half a pe- 
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FIG. 8: (color online). Close-packed structure of particles 
(large light circles) on a dual pyrochlore lattice. Thick (thin) 
lines represent tetrahedra from the lower (upper) kagome 
layer. Small dark circles denote vertices of one truncated 
tetrahedron. 

riod, which corresponds to the degeneracy A^'deg ~ 2^ . 
Once all chains are in phase, i.e., all particles occupy 
same vertex of the unit cell tetrahedron, we recover the 
initial structure constructed by filling hexagons in paral- 
lel kagome planes. 

To prove that the estimate ~ 1/12 gives the maxi- 
mal possible density of localized magnons we consider a 
polyhedron formed by 4 hexagons and 4 triangles, which 
has 12 vertices marked by small circles in Fig. |S1 This 
polyhedron is called truncated tetrahedron and belongs 
to the thirteen Archimedean solids. A special role of this 
polyhedron for our problem is determined by its geomet- 
rical structure: distance between two arbitrary vertices 
of a truncated tetrahedron is equal to d, \/3d, 2(i, or 
v^d. If one vertex of a truncated tetrahedron is occu- 
pied by a particle, then, according to the above exclusion 
rules all other vertices must be empty. Simple considera- 
tion yields that an A^-site pyrochlore lattice contains A/2 
truncated tetrahedra such that every lattice site is shared 
between six of them. Counting now becomes straight- 
forward: the number of hard-core particles in the effec- 
tive lattice-gas model cannot exceed one particle per six 
truncated tetrahedra (1/6) x {N/2), which yields 1/12 
as the upper bound on n. Above, we have constructed 
an explicit example of the particle arrangement with the 
density no — 1/12, therefore, this is indeed the highest 
density of the close-packed structure. 



B. Effective lattice-gas model 

Finite-temperature properties of the effective lattice- 
gas model have been studied by Monte Carlo simula- 
tions. In the low-density disordered regime z < 1 we 
use the same annealing protocol as in the 2D case and 
investigate cubic clusters with periodic boundary condi- 



tions and A^ = AL'^ sites, L = 6-12. At /x = (z = 1) the 
entropy obtained by numerical integration of MC data 
and the density of particles are equal to 

5/Ar = 0.1329(1) , n = 0.04718(1) , (12) 

which yields M — 1/2— n w 0.4528 for the magnetization. 

In order to study the high-density regime z > 1 we 
employ the exchange MC algorithm with 80 replicas dis- 
tributed in the range < Inz < 7. In contrast with 
the 2D case we still find significant equilibration problem 
for large lattices, which exhibit extremely long relaxation 
times t > 10^ MC steps despite of a high exchange rate 
^ 0.8 between replicas. Sufficient statistics has been ob- 
tained only for three system sizes: two cubic clusters 
with L = 4,6 and a tetragonal prism with 6 x 6 x 12 
unit cells (1728 sites). They were equilibrated for 5 x 10^ 
exchange MC steps and further 10^ MC steps were used 
for measurements. Evolution of the particle density with 
increasing z is presented in Fig. [5] for the two biggest 
systems. In the vicinity of Inz « 5.5 there is an abrupt 
jump of n to the value, which is very close to the density 
of the close-packed structure no = 1/12. 

We expect that the ordered state at Inz > 5.5 exhibits 
only a partial breaking of the translational symmetry. 
Similar to the studied 2D model, the half-filled chains re- 
main uncorrelated. This tendency should be even more 
pronounced for the pyrochlore lattice gas in view of an in- 
creased entropic effect from disordering determined by a 
larger A^dcg- As a result, the symmetry breaking includes 
only (i) selection of a chain direction and (ii) formation 
of \/3 x \/3 structure on a triangular lattice in the per- 
pendicular plane. In view of small available system sizes 
we have not performed investigation of the correspond- 
ing order parameters. Instead we have tried to clarify the 
nature of a possible phase transition into a close-packed 




FIG. 9: (color online). Density of a lattice gas on a pyrochlore 
lattice as a function of In 2: = fi/T. Arrow indicates the den- 
sity of a close-packed structure no = 1/12. Inset: probability 
distribution for the particle density obtained for A*" = 1728 
cluster at \nz = 5.2. 
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structure by collecting histograms for the particle den- 
sity. Results for N — 1728 lattice at Inz = 5.2 obtained 
with 10^ MC steps are presented in the inset of Fig. [9l 
A clear double peak structure with approximately equal 
weights in each peak points at a first-order transition. 
This should be, of course, confirmed by a similar study 
of larger clusters. We, therefore, put a tentative estimate 
In Zc = 5.2 ± 0.2 for the first-order transition into the 
magnon crystal state for a pyrochlore antiferromagnet. 

Since the original and the dual pyrochlore lattices have 
the same number of sites, the properties of a spin- 1/2 
pyrochlore antiferromagnet are obtained from the above 
results without any rescaling. In particular, the esti- 
mate for the entropy at the saturation field is given by 
Eq. (jl2|) . At zero temperature and H — Hs the magneti- 
zation jumps from M — 5/12 in the plateau region to the 
full saturation with M — 1/2. At finite temperatures the 
jump in the magnetization is much smaller AM « 0.015, 
see Fig. IHl Afterward, M{H) increases smoothly and 
asymptotically saturates. The first-order transition into 
the magnon crystal state takes place at 

Hc{T)^Hs-TlnZc^4:J-5.2T (13) 

for a spin- 1/2 model. 

IV. EFFECT OF EXTRA EXCHANGES 

In the previous two sections, we have found that the 
magnon crystal states in checkerboard and pyrochlore an- 
tiferromagnets exhibit massive degeneracy at T = 0. The 
degeneracy is not lifted at finite temperatures by fluc- 
tuations within the LM subspace. This leads to inter- 
esting types of broken symmetries in 2D and 3D cases. 
In principle, fluctuations to higher energy states outside 
the LM ensemble at T > can favor one or the other 
type of close-packed magnon structures. Such a scenario 
cannot be excluded on general grounds. The higher en- 
ergy states are, however, separated by finite gaps from 
LMs. Therefore, there must be a region near T ^ 0, 
H ^ Hs — 0, where such a selection is ineffective and 
only a partial breaking of translational symmetry takes 
place. Another mechanism of degeneracy lifting can be 
provided by the magnetoelastic coupling, though the 
corresponding study has considered only one possible 
pattern of LMs for each of the two lattices. In the follow- 
ing we investigate the third possibility: effect of further 
neighbor exchanges;^ which play significant role in some 
magnetic materials with pyrochlore lattices4iii^i^i^i^ 

We begin with the checkerboard antiferromagnet, 
considering an additional weak diagonal exchange of 
strength J', see Fig. [TId. The magnon crystal state be- 
comes completely unstable at a certain critical value of 
J'J^ We are interested in the subcritical regime | J'| 
0. Let us consider two localized magnons ^ on ad- 
jacent square plaquettes along the [11] direction: \ipi) 
and |</?i+2ai), see Fig. [5] and Sec. Ila for the notations. 
These LMs are connected by an additional exchange bond 



V = J'S2,i ■ S4^i+2ai- To measure energy shifts relative 
to the energy of the fully polarized state, we shall always 
subtract a constant from the exchange bond operators: 

V V - J'S^. The potential energy U of two LMs 
on adjacent voids is found by calculating the expectation 
value of V over the two-magnon state = \ipi(pi+2aLi) 
and correcting it by the self-energy of LMs without neigh- 
bors: 

U^{^P\V\^P)-2{^,\V\ip,)^J'/16. (14) 

This result applies to an arbitrary spin value S. For 
J' > presence of LMs on adjacent square voids is ener- 
getically unfavorable, because two spin flips with a flnite 
probability occupy the same J' bond. Degeneracy of the 
magnon crystal is, consequently, lifted in favor of an or- 
thorhombic structure: half-filled rows of LMs alternate 
in phase in order to minimize the contribution (jl4p be- 
tween the rows. For J' < localized magnons attract 
each other and the (2 x 2) tetragonal close-packed struc- 
ture is stabilized. 

The interaction energy between two LMs residing on 
hexagon voids of pyrochlore lattice is obtained in the 
same way as Eq. with the result U = J'/36. For 
a pyrochlore antiferromagnet three additional exchanges 
may be present :-" ' -5 i^j^g second-neighbor exchange J2 
for spin pairs at distance y/3d, the third-neighbor ex- 
change J31 for spin pairs at distance 2c? on opposite ver- 
tices of hexagons, and the third-neighbor exchange J32 
for spin pairs at distance 2d along chains, see Fig. [T^. 
(Details about the crystal structure of pyrochlore ox- 
ides can be found, e.g., in Ref. l38l .) Similar to the 
2D case the weak additional exchanges select relative 
shift between adjacent half-filled chains of LMs. If two 
chains are not shifted relative to each other, LMs occupy 
hexagons in parallel kagome planes, whereas if the shift 
is present, LMs form a nonplanar structure. Analysis 
of further neighbor exchange links between hexagons on 
a pyrochlore lattice is straightforward but cumbersome. 
Below we only summarize the conclusions. Once the two 
chains of LMs form a planar structure, one LM from the 
first chain is coupled to two hexagons in the second chain: 
to a hexagon in the same kagome plane by two J2 and 
two J32 bonds and to a hexagon in a parallel kagome 
plane by a J31 bond. Then, the interaction energy be- 
tween chains normalized per one LM of the first chain is 
Ui — (2J2 -I- J31 -t- 2J32)/36. In the nonplanar structure 
a given LM from the first chain is coupled to only one 
non-planar hexagon by three J2 bonds and one J32 bond, 
which yields U2 = (3J2 + •/32)/36. 

Comparing the interaction energies Ui and U2 for the 
two structures, we conclude that for ( J31 -I- J32 — J2) < 
localized magnons occupy hexagons in parallel kagome 
planes. Between the planes LMs follow the ABCABC... 
structure of close-packed hard spheres. If the above ex- 
pression changes sign, then LMs form a more complicated 
nonplanar structure. Amazingly enough, it is again frus- 
trated and degenerate: half-filled chains form a triangu- 
lar lattice, therefore, antiphase shifts would correspond 
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to an effective Ising antiferromagnet on a triangular lat- 
tice. The residual degeneracy should be finally lifted by 
zero-point fluctuations induced by the next-neighbor ex- 
changes. The corresponding energy scale is of the order 
of J'V J. 

V. SUMMARY 

Knowledge of the exact ground state and the single- 
particle excitation spectrum of frustrated quantum anti- 
ferromagnets above the saturation field Hs allows to de- 
velop a quantitative description of their thermodynamic 
properties in a finite range of fields H ^ Hg and tem- 
peratures T ^ 0. Localized nature of the lowest energy 
excitations leads to a mapping onto lattice gas models 
of hard-core classical particles with an appropriate ex- 
clusion principle. As a result, frustrated quantum anti- 
ferromagnets preserve a macroscopic degeneracy (finite 
entropy) at H — Hg and T = 0. The close-packed struc- 
ture of particles (localized magnons) corresponds to the 
magnon crystal state, which breaks only translational 
symmetry and preserves continuous rotations about the 
field direction. 

The main finding of our study is high degeneracy of the 
magnon crystal states in pyrochlore and checkerboard an- 
tifcrromagnets at zero temperature. This, however, does 



not exclude presence of finite-temperature phase transi- 
tions. The corresponding lattice gas models have been 
studied by the exchange Monte Carlo method. Specif- 
ically, we have found numerically the location Inzc — 
4.56(2) and two critical exponents /3 = 0.13(1) and 
f — 0.86(2) for a square-lattice gas with the second- 
neighbor exclusion, which describes a checkerboard anti- 
ferromagnet. Further numerical work is needed to clarify 
small differences with the recent independent MC study 
of the same modeli^ Monte Carlo simulations for the 
pyrochlore-lattice gas show a sharp jump in the parti- 
cle density, which suggests a first-order transition into a 
partially ordered crystal phase. 
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